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Introduction: 

The  purpose  of  this  project  is  to  develop  a  zebrafish  model  of  type  1  neurofibromatosis.  The  reason  that 
this  animal  model  may  be  useful  for  the  further  understanding  of  this  disease  is  that  genetic  and 
chemical  screens  will  be  possible  in  a  high  throughput  format  using  this  vertebrate  system.  In  addition, 
structure-function  analysis  will  be  possible  using  the  developing  zebrafish  model.  Prior  to  this  work,  the 
Nfl  genes  in  zebrafish  had  not  been  described  or  characterized.  Zebrafish  develop  rapidly  and  are 
transparent  during  development  allowing  for  easy  visualization  of  many  structures  including  the 
cardiovascular  system.  Genetic  and  chemical  screens  can  be  performed  in  high  throughput  modalities  at 
relatively  low  cost.  Emerging  technologies  allow  for  complex  manipulation  of  gene  expression.  This 
project  involves  the  collaboration  of  two  established  groups  to  coordinate  efforts  to  take  advantage  of 
zebrafish  as  a  model  system  for  NF1. 


Body  of  work: 

Identification  of  nfla  and  nflb 

We  utilized  a  bioinformatics  approach  to  identify  the  zebrafish  orthologues  of  human  NF1.  Analysis  of 
the  seventh  assembly  (Zv7)  of  the  zebrafish  genome  revealed  two  genes  that  are  highly  similar  to  NF1  at 
the  amino  acid  level  (90.4%  and  90.7%  respectively),  which  we  named  nfl  a  and  nflb.  These  genes  are 
highly  related  to  one  another  (87.4%  identical  and  93.7%  similar);  nfla  and  nflb  share  similar  genomic 
structures  and  each  contains  57  exons,  nfla  is  located  on  chromosome  15  (Figure  1A)  and  predicts  a 
311  kD  protein  composed  of  2755  amino  acids  while  nflb  is  located  on  chromosome  10  (Figure  IB)  and 
predicts  a  310  kD  protein  composed  of  2747  amino  acids. 

nfla  and  nflb  are  broadly  expressed  early  and  become  restricted  to  the  nervous  system  during 
later  embryogenesis. 

We  searched  the  public  database  for  zebrafish  Expressed  Sequence  Tag  (EST)  clones  using  the 
sequences  identified  from  Ensembl  and  found  many  ESTs  spanning  different  regions  of  the  nfla  and 
nflb  coding  sequences,  indicating  that  both  genes  are  expressed.  Among  them,  we  obtained  one  EST 
clone  for  each  gene,  in  order  to  evaluate  the  expression  of  nfla  and  nflb  during  development  (see 
Materials  and  Methods).  The  expression  of  nfla  and  nflb  throughout  early  blastula  stages  demonstrates 
the  maternal  contribution  of  these  genes  (Fig.  2A,  2F),  and  ubiquitous  expression  throughout  the  embryo 
continued  up  to  24  hours  post  fertilization  (hpf),  with  levels  being  highest  in  the  head  (Fig.  2B,  2G).  By 
2  days  post  fertilization  (dpf),  nfla  expression  remained  high  in  the  brain  (Fig.  2C-E),  as  demonstrated 
in  transverse  sections  along  the  diencephalon  and  tegmentum  (Fig.  2N),  unlike  the  expression  of  nflb , 
which  rapidly  decreased  with  time  (Fig.  2H-J).  Importantly,  both  nfla  and  nflb  continued  to  be 
expressed,  albeit  weakly,  throughout  the  spinal  cord  during  later  stages  (Fig.  2D-E  for  nfla;  2I-J  for 
nflb;  K  and  F  for  higher  magnifications),  nfla  expression  in  the  transverse  section  through  the 
embryonic  trunk  shows  expression  in  the  cytoplasm  of  both  neurons  and  glia  (arrowheads  in  Fig.  2M) 
and  is  clearly  evident  in  the  axonal  tracts  of  the  spinal  cord  (arrows  in  Fig.  2M).  nflb  expression  shows 
a  similar  pattern  of  expression  in  the  spinal  cord  but  is  much  weaker  (data  not  shown).  In  addition, 
strong  expression  of  nfla ,  but  not  nflb ,  is  detected  in  the  retina  during  later  stages  of  development  (Fig 
2D-E,  21- J).  Transverse  sections  show  that  nfla  is  enriched  in  the  ganglion  cell  layer  and  inner  nuclear 
layer  of  the  retina  (rgc  and  ini  in  Fig.  20),  and  nflb  is  weakly  expressed  in  the  gut  (Fig.  2P)  supporting 
the  possibility  that  these  genes  have  distinct  functional  roles  in  these  other  tissues.  However,  during 
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these  late  embryonic  stages,  the  principal  expression  domains  of  nfla  and  nflb  reside  in  the  central 
nervous  system  (CNS). 

The  number  of  OPCs  is  increased  upon  Nfl  knockdown 

In  order  to  test  the  role  of  nfla  and  nflb  in  zebrafish  oligodendrocyte  development,  we  used  the 
morpholino  knockdown  of  each  gene  in  the  olig2-EGFP  transgenic  zebrafish  line  bred  into  both  wild 
type  and  p53  -/-  backgrounds.  The  olig2-EGFP  transgenic  line  expresses  EGFP  in  OPCs  throughout 
oligodendrocyte  development.  The  olig2-EGFP  transgenic  line  in  the  p53  mutant  background  was  used 
to  circumvent  non-specific  effects  due  to  morpholino  toxicity.  To  knockdown  Nfl  function,  we 
designed  and  injected  nfl- specific  splice-blocking  morpholinos  (two  were  generated  against  nfla  and 
one  for  nflb )  into  the  olig2-EGFP  transgenic  zebrafish  lines.  The  efficiency  of  the  morpholinos  was 
confirmed  by  RT-PCR  at  3  dpf,  and  the  production  of  aberrant  bands  upon  morpholino  knockdown 
demonstrated  inappropriate  splicing  (arrowheads  in  Fig.  3 A;  morpholinos  nfla-e 7  and  nflb-e, 4  are 
shown).  The  fact  that  nfla-el  injection  led  to  several  aberrant  bands  suggested  that  there  were  multiple 
cryptic  donor  sites  within  the  exon  7  that  can  lead  to  splicing  of  nfla.  These  aberrant  bands  were 
sequenced  and  confirmed  to  contain  deletions  of  the  targeted  exons  or  insertions  of  additional  introns 
that  resulted  in  frame- shifts  in  coding  sequence  leading  to  early  transcript  termination  (data  not  shown). 

In  olig2-EGFP  transgenic  zebrafish,  a  subset  of  GFP+  OPCs  that  are  born  in  the  ventral  spinal  cord 
migrate  dorsally  and  differentiate  into  oligodendrocytes  (brackets  in  Fig.  3B,  3F,  3J),  while  a  second 
GFP+  OPC  subpopulation  remains  in  the  ventral  spinal  cord,  intermingled  with  olig2-EGFP  + 
motoneurons  and  interneurons.  When  nfla  and  nflb  were  knocked  down  using  morpholinos,  we 
consistently  found  an  increased  number  of  dorsally  positioned  OPCs  at  3  dpf  (Fig.  3B,  3F,  3J).  The  most 
dramatic  OPC  increase  relative  to  control  was  observed  when  both  the  nfla  and  nflb  genes  were 
knocked  down  together  in  the  p53  -/-  background  (Fig.  3J),  although  increases  were  also  observed  in 
individual  nfl  gene  knock-downs  (Fig.  3F,  3N).  These  results  indicate  that  both  nfl  genes  normally 
regulate  the  numbers  of  OPCs,  and  that  nfl  knockdown  results  in  an  increased  number  of  dorsally 
positioned  OPCs. 

This  phenotype  was  extended  to  the  evaluation  of  ventrally  positioned  OPCs  by  co-immunostaining 
with  an  anti-soxlO  antibody  that  labels  OPCs,  but  not  olig2-EGFP  +  motoneurons  and  interneurons  that 
are  located  in  the  spinal  cord  at  this  stage.  The  increased  number  of  olig2-EGFP  +;  soxlO-i-  OPCs  in 
nfla  and  lb  morphants  was  evident  in  transverse  sections  of  the  spinal  cord  (Fig.  3K-M,  compared  to 
3G-I).  The  mean  number  of  OPCs  was  determined  for  nfla+nflb  morphants  and  control  embryos  by 
counting  OPCs  in  immunostained  transverse  sections.  The  nfl  morphants  were  found  to  have  36%  more 
OPCs  than  control  p53  -/-  embryos  (9.0  vs.  6.6  respectively,  n=23  and  26,  p<0.0001;  Fig.  3N).  These 
findings  document  that  one  of  the  functions  of  nfl  is  to  negatively  regulate  the  number  of  OPCs  during 
development.  Although  the  OPCs  in  nfl  morphants  exhibit  signs  of  increased  motility  and  migration 
(see  below),  the  observed  increase  in  OPCs  is  not  due  to  an  aberrant  dorsal  migration  of  OPCs  that 
normally  would  remain  the  ventral  spinal  cord.  The  increase  in  OPC  numbers  was  observed  in  both 
migrating  and  non-migrating  OPC  populations  when  they  were  assessed  separately  (n=ll,  p<0.05). 
Also,  while  there  is  a  slightly  higher  increase  in  OPCs  in  the  p53  mutant  background,  presumably  due  to 
the  survival  of  cells  that  would  otherwise  succumb  to  morpholino  toxicity,  it  is  important  to  note  that 
significant  increases  in  OPCs  were  also  observed  in  the  wild  type  background  (Fig.  3N). 

nfl  loss  specifically  affects  OPCs  from  the  pMN  of  the  ventral  spinal  cord 

To  address  whether  the  loss  of  nfl  also  affected  olig2-EGFP- positive  motoneuron  derivatives 
forming  in  the  pMN  of  the  ventral  spinal  cord,  we  examined  the  expression  of  islet  1  which  is  expressed 
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in  primary/secondary  motor  neurons,  interneurons,  and  sensory  neurons  of  the  spinal  cord,  and  zn-5 
which  is  a  marker  for  developing  secondary  motoneurons.  Quantification  of  cell  numbers  using  anit-isl 
and  zn5  antibodies  failed  to  indicate  significant  differences  in  cells  expressing  either  of  these  markers  in 
nfla+b  knockdown  embryos  compared  to  control  s  (Fig.  4  F  vs  B;  N  vs  J;  Q  and  R)  .  In  addition, 
nfla+b  loss  did  not  affect  the  overall  size  or  shape  of  developing  secondary  motoneurons  (Fig  4  N  vs  J). 
These  results  suggest  that  nfl  knockdown  acts  specifically  on  the  oligodendrocyte  lineage,  while  the 
neuronal  lineages  from  the  pMN  domain,  such  as  motoneurons,  are  not  affected. 

nfl  loss  promotes  OPC  motility 

In  order  to  examine  the  cellular  behavior  of  OPCs  following  nfl  knockdown,  we  carried  out  live 
time-lapse  imaging  in  vivo.  We  monitored  the  olig2-EGFP  transgenic  animals  for  12  hours,  beginning 
at  60  hpf  when  OPCs  are  actively  migrating  away  from  the  ventral  spinal  cord.  Compared  to  OPCs  in 
uninjected  control  p53  -/-  animals,  more  OPCs  in  the  nfla+  //^-knockdown  p53-/-  embryo  were  observed 
to  enter  the  dorsal  spinal  cord  (Fig.  5A,  B),  consistent  with  our  observations  in  fixed  sections  (Fig.  3). 
The  time-lapse  study  also  revealed  that  the  motility  of  nfla+lb- morphant  OPCs  is  enhanced  (Fig.  5E). 
When  the  velocity  of  individual  migrating  OPCs  was  measured  based  on  the  total  distance  traversed 
over  time  (|im/5  minute),  a  36%  increase  was  observed  upon  nfl  knockdown  (0.2686  vs.  0.3667  pm/5 
min;  p<0.005;  Fig.  5E).  n(la+  //^-knockdown  OPCs  migrate  further  on  average  in  the  dorsal  direction 
than  control  OPCs.  They  also  show  greater  tangential  movement,  presumably  due  to  their  increased 
motility,  but  there  is  no  net  migration  in  either  the  rostral  or  caudal  directions  as  shown  by  the  individual 
migratory  traces  (Fig.  5  C  and  D).  Therefore,  nfl  loss  leads  to  an  enhancement  of  motility  in  OPCs,  a 
phenotype  that  has  not  been  previously  described. 

nfl  loss  does  not  activate  Ras  signaling  pathways  in  OPCs 

The  best-known  function  of  Nfl  is  to  act  as  a  Ras-GAP  to  down-regulate  Ras  signaling  through  its 
downstream  pathways.  To  investigate  the  effect  of  nfl  loss  on  these  pathways  we  evaluated  the  status  of 
1)  phosphorylated  ERK1/2,  a  downstream  target  of  the  Ras/Raf  pathway  and  2)  phosphorylated 
ribosomal  protein  S6,  a  downstream  target  of  the  PI3  Kinase  pathway.  At  3  dpf,  when  the  increased 
OPC  numbers  were  evident,  and  more  cells  in  transverse  sections  of  the  spinal  cord  nfl a+ lb- 
knockdown  embryos  were  labeled  with  a  phospho-ERKl/2  antibody  than  in  controls,  indicating  more 
cells  with  ERK  activation.  This  is  consistent  with  Nfl’s  role  as  a  suppressor  of  signaling  through  the 
RAS/MAPK  pathway  (Fig.  6).  However,  we  were  surprised  to  find  that  the  phospho-ERKl /2-positive 
cells  in  the  dorsal  spinal  cord  did  not  correspond  to  GFP-positive  OPCs  in  either  control  or  nfl- 
knockdown  conditions.  Thus,  although  the  RAS/MAPK  pathway  was  significantly  activated  in 
surrounding  neurons,  there  was  not  a  detectable  increase  in  activity  in  the  OPCs  themselves.  We  also 
examined  Akt/mTOR  pathway  activation  and  found  that  while  phospho-S6-positive  cells  were  detected, 
nfl  loss  did  not  cause  significant  changes  in  the  numbers  of  phospho-S6-positive  cells  in  the  dorsal 
spinal  cord,  and  GFP-positive  OPCs  labeled  weakly  if  at  all . 

Morpholino  knockdown  of  nfla  and  nflb  results  in  cardiovascular  defects 

nfl  a  and  nflb  morphant  embryos  also  displayed  gross  abnormalities  of  cardiovascular 
development  appreciable  to  the  blinded  observer  by  48  hpf  (Fig.  7).  For  example,  2  ng  of  an  nfla 
translation-blocking  MO  produced  cardiac  defects  in  75%  of  injected  embryos  at  48  hpf,  while  only  5% 
of  embryos  injected  with  a  control  MO  (containing  5  altered  nucleotides)  showed  defects  when  assessed 
by  a  blinded  observer.  Morphant  embryos  frequently  displayed  a  “to-and-fro”  movement  of  blood 
within  the  heart.  That  is,  blood  was  seen  to  move  back  and  forth  from  atrium  to  ventricle,  suggesting  a 
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malfunctioning  atrio- ventricular  valve.  In  addition,  we  observed  pooling  of  blood  in  the  common 
cardinal  vein  and  a  paucity  of  blood  flow  along  the  dorsal  aorta  and  posterior  cardinal  vein.  Valvular 
insufficiency  and  reduced  blood  flow  were  not  appreciated  in  control  morphants  or  wild  type  embryos. 
Overall  development  of  the  embryos  was  relatively  preserved  through  the  first  3  days  of  development 
despite  these  cardiac  defects.  Histological  analysis  revealed  a  thinned  ventricular  myocardium  and  large 
pericardial  effusions  in  MO-treated  embryos.  Gross  morphological  analysis  also  showed  an  increased 
incidence  of  pericardial  effusions  at  48  hpf,  reflecting  cardiac  dysfunction,  in  nfla  and  nflb  morhpants 
when  compared  with  controls  (Fig.  7).  Non-specific  toxicity  arising  from  MO  exposure  as  a  cause  of 
the  observed  cardiovascular  defects  was  unlikely  as  unrelated  control  or  scrambled  MOs  failed  to 
produce  similar  levels  of  abnormalities,  defects  were  observed  even  at  low  doses  of  specific  MOs,  and 
similar  defects  were  observed  with  several  unrelated  but  specific  MOs  directed  against  nfla  and  nflb. 

In  addition,  injection  of  specific  MOs  in  p53  mutant  embryos  also  produced  similar  cardiovascular 
defects,  and  off-target  effects  due  to  MO  exposure  is  known  to  be  at  least  partially  mediated  through  p53 
activation.  Defects  in  cardiac  valve  morphogenesis  and  a  thinning  of  the  ventricular  myocardium  are 
also  seen  in  Nfl -deficient  murine  embryos. 

We  performed  knockdown  experiments  using  zebrafish  embryos  in  which  endothelial  cells  are 
marked  by  expression  of  a  cytoplasmic  enhanced  green  fluorescent  protein  (EGFP)  in  order  to  allow  for 
a  more  detailed  analysis  of  vascular  development.  Dramatic  abnormalities  of  vascular  patterning  in  the 
intersomitic  vessels  of  morphant  embryos  were  present  at  48  and  72  hpf  (Fig.  8).  In  nfla  MO  treated 
embryos,  the  leading  edge  of  the  sprouting  vessels  displayed  claw-like  projections  at  48  hpf  (Fig.  8C) 
and  failed  to  pattern  normally  such  that  the  dorsal  longitudinal  anastomotic  vessel  (DLAV)  did  not  form 
or  developed  only  in  a  rudimentary  fashion  (Fig.  8F).  This  occurred  in  embryos  that  were  otherwise 
normal  in  overall  size  and  maturity.  Vascular  pattering  defects  did  not  appear  to  correlate  directly  with 
cardiac  defects,  as  we  observed  embryos  with  vascular  abnormalities  that  did  not  display  pericardial 
effusion  or  valvular  insufficiency  as  assessed  by  a  “to-and-fro”  movement  of  blood  within  the  heart. 
Blood  flow  within  the  dorsal  aorta  appeared  intact  in  these  embryos. 

At  earlier  time-points,  we  assessed  vascular  defects  caused  by  nfla  and  nflb  knockdown  in 
zebrafish  embryos  expressing  a  nuclear-localized  EGFP  in  endothelial  cells.  Morphant  embryos 
displayed  a  complete  (Fig.  9A4)  or  partial  absence  of  intersomitic  vessels  emanating  from  the  dorsal 
aorta  at  24  hpf  when  compared  with  stage-matched  controls  (Fig.  9A3).  Overall  morphology  of 
morphant  and  control  embryos  appeared  equivalent  (Fig.  9A1-2)  ruling  out  non-specific  developmental 
delay.  These  defects  were  apparent  following  MO-mediated  knockdown  of  nfla  or  nflb  while 
knockdown  of  both  together  had  an  additive  effect  (Fig.  9B).  Again,  similar  defects  were  observed  with 
several  unrelated  but  specific  MOs  directed  against  nfla  and  nflb  (Fig.  10).  Our  analysis  of  morphant 
embryos  at  24  hpf  also  revealed  a  caudal  vessel  defect.  Morphant  embryos  displayed  a  cystic  expansion 
in  the  region  of  the  caudal  artery  and  vein  (Fig.  10B2)  when  compared  with  controls  (Fig.  10B1). 
Identity  of  the  expanded  tissue  as  vascular  was  confirmed  by  expression  of  EGFP  (Fig.  10B3-4)  and  the 
presence  of  red  blood  cells  in  the  expanded  region.  This  defect  was  present  following  knockdown  of 
either  nfla  or  nflb  (Fig.  10C). 

Additional  confirmation  of  the  role  of  nfla  in  vascular  development  derives  from  studies  using  a 
genetic  background  sensitized  to  vascular  insult.  Flt4  is  the  zebrafish  orthologue  of  VEGF  receptor-3, 
and  previous  studies  have  employed  MOs  directed  against  flt4  to  investigate  genetic  interactions  during 
zebrafish  artery  development.  Additionally,  flt4  morphant  zebrafish  embryos  display  variable  defects  in 
segmental  artery  formation  reminiscent  of  those  identified  in  our  nfl  morphants.  As  flt4  expression  is 
restricted  to  endothelial  cells,  combined  nfla  and  flt4  knockdown  would  be  expected  to  reveal  genetic 
interactions  that  are  preferential  to  the  vasculature.  Endothelial-GFP  expressing  embryos  were  injected 
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with  /7/4  MO  alone  and  in  combination  with  a  MO  directed  against  nfla,  and  morphant  embryos  were 
assessed  at  24  and  48  hpf.  At  the  low  doses  of  MOs  used,  85%  of  flt4hif la  compound  morphants 
displayed  abnormal  vascular  shunts  at  48  hpf  as  compared  to  only  4-8%  of  individual  flt4  or  nfla 
morphants  (Figure  9C).  This  defect  was  not  apparent  in  control  morphants. 

Zebrafish  nfla  and  nflb  mutant  lines 

A  ciritical  part  of  our  proposal  involves  the  generation  of  stable  zebrafish  lines  carrying 
mutations  in  nfla  and  nflb.  These  animals  will  verify  our  morpholino  studies  and  be  able  to  extend  our 
analysis  beyond  transient  assays.  Furthermore,  planned  future  studies  include  chemical  and  genetic 
screen  strategies  based  on  these  animals  that  can  contribute  to  therapeutic  development.  Through 
collaborative  efforts  we  have  used  a  reverse  genetics  technology  called  TILLING  (Target  Induced  Local 
Lesions  In  Genomes)  and  found  several  nfla  missense  mutant  alleles  and  one  nonsense  mutant  line  in 
the  nfla  gene  (Fig.  1 1A  and  B).  A  second  mutagenesis  strategy  based  on  retroviral  insertion  identified  a 
single  mutant  line  affecting  the  nflb  locus  in  the  intron  between  exon  1  and  2  (Fig.  11B).  These  nfla 
and  nflb  mutant  alleles  have  recently  been  recovered  in  our  laboratory  and  being  currently  being  tested 
for  potential  defects  in  oligodendrocyte  and  cardiovascular  development. 

Most  recently,  a  new  tool  was  developed  allowing  targeted  gene-specific  mutagenesis  using 
Zinger  Finger  Nuclease  (ZFN)  technology  in  the  zebrafish  research  field  ’ "  .  Briefly,  ZFN  is  an 
artificial  nuclease  designed  to  bind  specific  nucleotide  sequences  in  a  zinc  finger-based  domain  and 
induced  to  cleave  double  strand  DNA  near  the  target  sequence.  Non-homologous  end-joining  repair 
results  in  small  deletions  or  insertions  in  the  gene  of  interest.  We  designed  several  ZFNs  that  can  induce 
mutations  in  the  specific  sequences  of  the  nfla  and  nflb  genes  and  obtained  three  frame- shift  mutations 
(. nflam ,  nfla 053  and  nfla D59)  and  one  in-frame  insertion  mutation  ( nfla+33ID& )  in  exon  25  of  nfla  and  an 
in-frame  deletion  mutation  ( nflbD 39)  in  exon  17  of  nflb.  Because  the  GRD  domain  in  nfla  and  nflb 
spans  exon  26  and  exon  27,  we  expect  that  three  frame- shift  mutant  alleles  for  nfla  should  lack  GAP 
activity.  Our  initial  preliminary  study  examining  oligodendrocyte  development  in  one  of  these  lines  is 
very  promising.  We  have  bred  the  mutant  carrying  an  8bp  deletion  in  exon  25  of  nfla  ( nfla  )  to  the 
olig2:gfp  transgenic  line.  The  homozygous  nflam  mutant  embryos  appear  morphologically  normal  and 
survive  through  two  weeks  of  age  indicating  that  they  may  be  viable  through  adulthood.  Importantly, 
we  observed  increased  number  of  GFP+  OLCs  in  nflam-;Tg(olig2:gfp)  mutant  embryos  when 
compared  to  wild-type  or  heterozygous  nflam  embryos  (Fig.  12).  We  have  also  assayed  oligdendrocyte 
development  in  the  nflb  in-frame  mutant,  however  no  defects  or  differences  in  oligodendrocyte 
development  have  been  detected  thus  far  (data  not  shown).  We  are  also  just  beginning  to  use  a  modified 
ZFN  strategy  that  is  reported  to  have  improved  efficacy,  and  we  are  employing  this  approach  to  further 
insure  our  obtaining  stable  mutants  in  the  nflb  gene. 

In  summary 

•  The  number  of  OPCs  in  the  spinal  cord  are  increased  upon  transient  knockdown  of  nfla  and  nflb 

•  Phospho  MAP  kinase  pathway  may  not  be  responsible  for  the  OPC  increase  or  be  required  non-cell 
autonomously 

•  Cardiac  and  vascular  defects  in  nfla  and  nflb  knockdown  fish 

•  Several  putative  mutants  for  nfla  and  nflb  have  been  identified  and  validated 

Materials  and  Methods 

1.  Immunohistochemistry  Fixed  embryos  were  embedded  in  sucrose/agar  mix  and 
cryosectioned  using  cryostat  machine  with  14pm  thick.  Sections  were  incubated  with  the  primary 
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antibodies  (soxlO,  1:5,000;  GFP,  1:500,  pMAPK,  1:200,  pH3,  1:200)  and  Alexa  568  or  633  secondary 
antibodies. 

2.  Whole  mount  in  situ  hybridization  Fixed  embryos  were  permeabilized  with  proteinase  K 
treatment  (lOug/ml)  with  various  incubation  time  depending  on  the  stages  of  embryos,  which  then  were 
incubated  with  DIG-labeled  antisense  RNA  probes  and  visualized  using  alkaline  phosphatase- 
conjugated  anti-DIG  antibody  in  color  substrates. 

3.  RT-PCR  cDNA  for  RT-PCR  was  made  by  reverse  transcribing  RNA  prepared  from 
morpholino  (MO)  injected  embryos  using  Trizol,  and  PCR  was  performed  to  detect  potential  aberrant 
splicing.  The  triplicates  of  the  same  cDNA  were  used  for  the  quantitative  PCR  to  measure  the 
knockdown  by  nflb-e4  with  beta-actin  as  baseline  control. 

4.  Confocal  imaging  Zeiss  LSM  510  META  confocal  microscope  at  Harvard  Neurodiscovery 
imaging  program  was  used  for  capturing  confocal  images  of  olig2-EGFP  and  other  antibody  staining  in 
zebrafish  spinal  cord  with  l-3pm  optical  sections. 

5.  Counting  OPCs  and  statistical  analysis  Three  sections  of  the  spinal  cord  at  the 
approximately  similar  anterior-posterior  level  from  individual  embryos  (eight  to  ten  embryos)  were 
chosen  for  counting  all  of  the  GFP+/soxl0+  cells  within  the  section.  Individual  t-test  were  performed  to 
test  the  significant  differences  in  the  number  of  GFP+/soxl0+  cells  in  the  spinal  cord  using  Prism 
software. 

6.  TUNEL  Chemicon  ApopTag  Red  In  Situ  Apoptosis  Detection  Kit  (S7165)  was  used  for 
TUNEL  staining.  The  protocol  for  cryosections  in  the  munual  was  followed  with  minor  modifications. 
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Fig.  1.  Zebrafish  have  two  orthologues  of  human  NF1. 

(A.  B)  Genomic  and  mRNA  structures  of  the  two  orthologous  zebrafish  genes  corresponding  to  human 
NF1.  (C)  Phylogenetic  comparison  of  zebrafish,  human,  mouse,  and  Drosophila  neurofibromin.  (D) 
Analysis  of  syntenic  relationships  between  human  chromosome  17  ( NF1 )  and  zebrafish  chromosomes 
15  (nfla),  and  10  (nflb).  Both  genes  are  flanked  by  upstream  genes  encoding  WD  repeat  and  SOCS 
box-containing  1  (WSB1),  Kinase  suppressor  of  Ras  1  (KSR1),  and  Galectin-9  (LGALS9)  and 
downstream  genes  encoding  A  kinase  anchor  protein  1  (AKAP1)  and  RNA-binding  protein  Musashi 
homolog  2  (MSI2).  Relative  genomic  positions  are  to  scale  as  indicated. 
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Fig.  2.  The  expression  of  nfla  and  nflb  during  embryogenesis  by  whole  mount  in  situ 
hybridization. 

nfla  expression:  panels  A-E,  K,  M,  N,  O.  nflb  expression:  panels  F-J,  P.  A,  F,  one  cell  stage;  B,  G, 
24hpf;  C,  H,  48hpf;  D,  I,  3dpf;  E,  J,  5dpf.  (A-E)  nfla  RNA  is  maternally  supplied  and  expressed 
strongly  throughout  the  embryo.  Expression  remains  strong  in  the  brain  as  development  proceeds.  (F- 
J)  nflb  RNA  is  also  maternally  provided  and  expressed  broadly  at  a  high  level  during  early 
embryogenesis,  but  its  expression  drops  to  low  levels  in  the  whole  embryo  after  2dpf.  (K,  L)  At  3dpf, 
both  nfla  and  nflb  expression  is  weaker  and  becomes  restricted  to  the  spinal  cord.  (M)  Transverse 
section  of  the  spinal  cord  hybridized  with  the  nfl  a  probe  at  3dpf .  Expression  is  cytoplasmic  in  neurons 
and  glial  cells  (arrowheads),  and  axon  tracts  (arrows).  (N)  Expression  of  nfla  in  a  transverse  section 
through  a  3dpf-old  embryo  at  the  level  of  the  optic  chiasm  (oc),  showing  broad  expression  in  the 
diencephalon  (dc)  and  the  tegmentum  (tg).  (O)  Transverse  section  of  a  5dpf-old  embryo  retina 

hybridized  with  nfla,  showing  strong  expression  in  the  retinal  ganglion  cells  (rgc)  and  inner  nuclear 
layer  (ini),  onl,  outer  nuclear  layer.  (P)  nfla  and  nflb  are  also  expressed  outside  the  CNS.  nflb 
expression  in  a  cross  section  through  the  gut  at  5dpf  is  shown,  fb,  forebrain;  mb,  midbrain;  hb, 
hindbrain;  sc,  spinal  cord;  nc,  notochord.  Scale  bars  in  A-E  =  100pm;  K,  L,  N,  O  =  50p.m;  M,  P  =  25p.m. 
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Fig.  3.  The  knockdown  of  nfla  and  nflb  leads  to  increased  numbers  of  OPCs. 

(A)  RT-PCR  assays  following  the  injection  of  nfla-el  and  nflb-e 4  morpholinos  individually  or  in 
combination  into  olig2-EGFP  or  olig2-EGFP\p53-/-  embryos  at  3dpf.  The  table  above  the  picture 
denotes  the  injected  morpholinos.  Morpholino  injections  induce  aberrant  splicing  resulting  in  the 
production  of  abnormal  bands  indicated  by  arrowheads  while  the  arrows  indicate  wild  type  bands.  (B-E) 
Uninjected  olig2-EGFP  transgenic  zebrafish.  (F-I)  Control  morpholino  injected  olig2-EGFP;  p53-/- 
embryo.  (J-M)  olig2-EGFP\  p53-/-  injected  with  nfla  and  nflb  morpholinos.  (B,  F,  J)  Projected 
confocal  images  show  a  lateral  view  of  the  spinal  cord  in  live  transgenic  animals  at  3dpf.  The  number  of 
OPCs  is  increased  upon  morpholino  injection  (J),  compared  to  uninjected  or  control  morphants  (B,  F). 
(C-E,  G-I,  K-M)  Transverse  sections  through  the  spinal  cord  of  embryos  at  3dpf,  labeling  GFP  (green) 
and  soxlO  (red),  to  specifically  visualize  OPCs  (arrows).  The  number  of  GFP+/soxlO+  OPCs  in 
morphant  embryos  is  increased  compared  to  controls,  consistent  with  the  lateral  views  (B,  F,  J),  and 
confirms  that  OPCs  in  both  the  dorsal  (migrating)  and  ventral  (nonmigrating)  spinal  cord,  increase  in 
numbers.  (N)  Statistical  analysis  of  averaged  numbers  of  GFP+/soxlO+  cells  per  section  in  nfla  (la)  or 
nflb  (lb)  alone  or  in  combination  at  3pdf.  Error  bar  indicates  SEM,  asterisks  indicate  statistical 
significance. 
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Fig.  4.  Spinal  cord  neurons  are  not  affected  by  Nfl  loss. 

(A-D,  I-L)  olig2-EGFP\p53-/-  embryos  injected  with  the  control  morpholino.  (E-H,  M-P)  An  olig2- 
EGFP\  p53-/-  embryo  injected  with  nfl  a  and  nflb  morpholinos.  (A-H)  Projected  confocal  images  of 
transverse  sections  labeled  with  the  anti -isl  antibody  show  subsets  of  primary/secondary  motor  neurons 
(bracket  1),  interneurons  (bracket  2),  and  sensory  neurons  (asterisk  with  a  question  mark).  The  number 
and  morphology  of  anti-/.v/-positivc  cells  appears  to  be  unaffected  upon  nfl  knockdown.  Green,  olig2- 
GFP;  red,  anti -isl;  blue,  DAPI.  Arrows  in  A,  C,  E,  G,  dorsal  OPCs.  (IP)  Fluorescence  images  of 
transverse  sections  labeled  with  the  zn5  antibody  showing  secondary  motor  neurons.  The  number  and 
morphology  of  zn-5-positive  cells  appears  to  be  unaffected  upon  nfl  knockdown.  Statistic  analyses  of 
the  numbers  of  anti-isll  (Q)  or  zn5  (R)  show  no  significant  differences  (from  15  sections).  Green,  olig2- 
GFP;  red,  zn5;  blue,  DAPI.  Arrows  in  I,  K,  M,  O  indicate  dorsal  OPCs;  asterisks  in  J,  K,  N,  O  indicate 
the  dorsal  lateral  fasciculus;  brackets  in  J,  K,  N,  O  indicate  secondary  motoneurons.  Scale  bar  =  20pm. 
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Fig.  5.  nfl  loss  enhances  the  motility  of  OPCs 

(A,  C)  An  uninjected  olig2-EGFP;  p53-/-  control  embryo;  (B,  D)  An  olig2-EGFP\  p53-/-  embryo 
injected  with  nfla  and  nflb  morpholinos.  (A,  B)  Montages  of  12  hour  time  lapse  images  from  60hpf  to 
72hpf,  showing  the  movements  of  OPCs  every  hour.  Original  images  were  taken  every  five  minutes.  For 
the  complete  movies,  refer  to  Supplementary  video  1.  Compared  to  the  uninjected  control,  OPCs  with 
nfl  loss  generally  migrate  at  a  higher  velocity.  (C,  D)  Cell  migration  paths  are  shown  for  OPCs  that 
exhibited  the  highest  calculated  velocities.  In  either  condition,  their  migration  patterns  appeared 
stochastic.  Arrowheads  indicate  the  end  point  of  the  cell’s  migration.  Top,  dorsal;  left,  rostral.  (E) 
Graph  showing  the  migration  velocity  of  OPCs  with  or  without  Nfl  loss.  Each  triangle  represents  the 
velocities  of  ten  individual  OPC  cells  demonstrating  the  highest  motility.  n=10  for  uninjected  control 
embryos,  n=12  for  nfl a+ 1  b  morphant  embryos.  Scale  bars  in  A,  B  =  25pm.  Scale  bars  in  C,  D,  X  axis  = 
10pm,  Y  axis  =  5  pm. 
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Fig.  6.  Nfl  loss  activates  the  RAS/MAPK  pathway  in  the  spinal  cord,  but  not  in  OPCs. 

Projected  confocal  images  of  transverse  spinal  cord  sections  labeled  with  GFP-(green)  and  phopho- 
ERKl/2-specific  (red)  antibodies.  (A-D)  An  uninjected  olig2-EGFP\p53-/-  control  embryo.  (E-H) 
o  I  ig2  -  E  G  FP ;  p53  -/-  embryo  injected  with  nfl  a  and  nflb  morpholinos.  Both  the  intensity  of  phospho- 
Erkl/2  staining  and  the  number  of  phospho-Erkl/2-positive  cells  are  increased  upon  nfl  loss;  however, 
migrating  GFP+  OPC  cells  in  the  dorsal  spinal  cord  do  not  exhibit  phospho-ERKl /2-labeling  (arrows). 
Blue  DAPI  staining  indicates  cell  nuclei  (D,  H). 
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Fig.  7.  MO  knockdown  of  nfla  results  in  gross  defects  of  cardiovascular  development. 

Quantification  of  cardiovascular  defects  in  72  hpf  nfla  ATG  MO-  and  nfla  ATG  5MP  MO-injected 
(2ng)  embryos. 
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Fig.  8.  MO  knockdown  of  nf la  results  in  vascular  patterning  defects  at  48  and  72  hpf. 

(A-C)  At  48  hpf,  nfla  ATG  MO  treated  Tg(fli:egfp)yl  (endothelial- specific  GFP  transgenic)  zebrafish 
embryos  display  gross  defects  in  vascular  development  as  compared  with  control  MO  treated  or 
uninjected  samples.  Morphant  embryos  (C)  display  abnormal  projections  at  the  leading  edge  of  the 
developing  intersomitic  vessels  and  fail  to  develop  the  dorsal  longitudinal  anastomotic  vessel  (DLAV) 
present  in  both  controls  MO  treated  (B)  and  uninjected  (A)  samples  (arrows).  (D-F)  At  72  hpf,  nfla 
ATG  morphant  embryos  display  only  rudimentary  DLAVs  (arrow)  and  disorganization  of  the 
vasculature  (F)  when  compared  with  control  MO  treated  (E)  or  uninjected  (D)  embryos.  Scale  bars: 

25  pm. 
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Fig.  9.  MO  knockdown  of  nfla,  nflb,  or  both  result  in  vascular  defects  at  24  hpf. 

(A,  B)  Analysis  and  quantification  of  vascular  defects  at  24  hpf  in  uninjected  and  morphant 
Tg(fli:negfp)y7  (endothelial- specific  nuclear  GFP  transgenic)  zebrafish  embryos.  Control  MO-treated 
(Al)  and  combined  nflci/ nflb  MO-treated  (A2)  zebrafish  embryos  appear  similar  by  gross 
morphological  analysis  at  24  hpf  (scale  bars:  500pm).  Development  of  intersomitic  vessels  is  deficient 
at  24  hpf  in  nfla/ nflb  combined  morphants  (A4)  when  compared  to  control  MO-treated  embryos  (A3) 
(scale  bars:  25pm).  (B)  Intersomitic  vessel  formation  between  somites  17-30  at  24  hpf  was  scored  as 
absent  (red),  intermediate  (grey)  or  normal  ( nfla  ATG  5MP  MO  refers  to  a  control  MO  with  5  mispairs 
when  compared  to  the  nfla  ATG  MO).  (C)  MO-mediated  knockdown  of flt4,  providing  a  sensitized 
background  for  the  detection  of  vascular  defects,  was  combined  with  nfla  ATG  MO  knockdown.  85% 
of  combined  flt4/nfl a  MO-treated  embryos  display  abnormal  vasculature  as  compared  with  4-8%  of 
individual  flt4  or  nfla  MO-treated  embryos. 
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Fig.  10.  MO  knockdown  of  nfla,  nflb,  or  both  together  result  in  vascular  defects  at  24  at  48  hpf. 

(A)  Quantification  of  vasculature  at  24  hpf  in  uninjected  and  morphant  Tg(fli:negfp)y7  (endothelial- 
specific  nuclear  GFP  transgenic)  zebrafish  embryos,  nfla  morphant  embryos  were  qualitatively 
assessed  for  presence,  absence  (red),  or  an  intermediate  phenotype  (grey)  with  regard  to  the  developing 
trunk  intersomitic  vessels  between  somites  17-30  at  24  hpf.  (B 1-B4)  Treatment  of  Tg(fli:negfp)y7 
zebrafish  embryos  with  an  nfla  ATG  MO  also  leads  to  a  cystic  expansion  in  the  region  of  the  caudal 
artery  and  vein  (B2)  not  present  in  uninjected  embryos  (B 1 )  (scale  bars  represent  100pm).  The 
expanded  tissue  was  confirmed  to  be  vascular  by  GFP  expression  (B3,  B4)  (scale  bars  represent 
100pm).  (C)  This  phenotype  is  present  in  nfla  ATG,  nflb  ATG,  and  nfla/nflb  combined  ATG  MOs  as 
quantified. 
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Fig.  11.  Searching  for  zebrafish  mutants  for  nfla  and  nflb 

(A)  Schematic  of  TILLING  procedure.  (B)  A  list  of  the  mutations  for  nfla  and  nflb  identified  from 
TILLING  thus  far. 
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Fig.  12.  OPC  phenotype  is  recapitulated  in  nfla  mutant  zebrafish  line  generated  using  engineered 
Zinc  Finger  Nucleases. 

Increased  number  of  GFP  OLCs  in  nfla  ;Tg(olig2:gfp)  mutant  embryos  (B)  when  compared  to 
wild-type  (A)  or  heterozygous  nflam  embryos. 
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Key  research  accomplishments: 

•  Characterization  of  zebrafish  orthologs:  nfla  and  nflb. 

•  Analysis  of  gene  expression  of  nfla  and  nflb. 

•  Analysis  of  morpholino  knockdown  phenotypes  for  nfla  and  nflb. 

•  Identification  of  fish  with  mutations  in  nfla  and  nflb  via  TILLING  and  through  the  use  of 
engineered  Zinc  Finger  Nucleases. 

•  Obtained  and  began  breeding  fish  with  insertional  mutation  in  nflb. 

Reportable  outcomes: 

Two  manuscripts  describing  the  results  obtained  so  far  have  been  submitted  for  publication  and  are 
presently  being  revised  to  address  reviewers’  concerns. 

Conclusion: 

In  this  second  year  of  study,  we  have  furthered  our  advances  towards  meeting  the  goals  set  forth  in  the 
funded  proposal.  The  target  genes  have  been  clearly  identified  and  phenotypes  due  to  morpholino 
knockdown  have  been  characterized  that  are  relevant  to  the  understanding  of  human  disease.  While 
morpholino  knockdown  is  informative,  future  experiments  will  be  aided  by  the  characterization  of  stable 
mutations  within  these  genes  in  zebrafish.  Significant  strides  have  been  made  towards  obtaining  these 
mutant  fish  lines.  Several  different  approaches  have  been  undertaken  which  may  result  in  an  allelic 
series  of  mutant  fish  appropriate  for  further  analysis  throughout  the  scientific  community.  When 
characterized,  these  reagents  will  be  freely  shared  with  our  colleagues.  These  reagents  should  be 
excellent  candidates  for  use  in  screening  small  molecules  and  for  use  in  performing  genetic  screens  to 
identify  modulators  of  this  important  disease. 
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